Three-dimensional imaging of colloidal glasses under steady shear 
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Using fast confocal microscopy we image the three-dimensional dynamics of particles in a yielded 
hard-sphere colloidal glass under steady shear. The structural relaxation, observed in regions with 
uniform shear, is nearly isotropic but is distinctly different from that of quiescent metastable colloidal 
fluids. The inverse relaxation time t^ 1 and diffusion constant D, as functions of the local shear rate 
7, show marked shear thinning with r" 1 oc D oc -y 0,8 over more than two decades in 7. In contrast, 
the global rheology of the system displays Herschel-Bulkley behavior. We discuss the possible role 
of large scale shear localization and other mechanisms in generating this difference. 

PACS numbers: 83.50.Ax, 83.60.-a, 83.80.Hj, 83.85.Ei 
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Glassy materials are ubiquitous in nature and in indus- 
try; examples range from molecular and metallic glasses 
[l|, 0] to soft glasses like colloidal suspensions, emulsions 
and foams [1, 0). Of special importance, both funda- 
mentally and practically, is their rheological behavior. 
Glasses have liquid-like microstructure, but solid-like me- 
chanical behavior. At low applied stress, they have finite 
shear moduli, but at sufficient stresses they yield and 
display highly nonlinear flow behavior. 

Among the many open issues in nonlinear glassy rhe- 
ology, steady shear stands as the simplest example,yet 
it is far from being fully understood. Theories 0,11, H0] 
have invoked various mechanisms for shear-induced re- 
laxation of initially arrested structures, predicting a va- 
riety of constitutive relations. Significantly, all these the- 
ories assume globally uniform shear. Simulations, so far 
the main tool to check directly the relation between mi- 
croscopic behavior and macroscopic flow, reveal spatially 
heterogeneous relaxation 0, Q and strong shear thinning 
d, Experiments are just starting to address micro- 
scopic dynamics under shear, but have been limited to 
coarse-grained data, two dimensional (2D ) or interrupted 




.HUGS. Moreover, 
as well as bound- 



flows, or ordering phenomena 
experiments imaging global flow 
ary driven simulations of Lennard- Jones (LJ) glasses [16| 
show that (soft) glasses often exhibit shear localization, 
which can not be described by simple constitutive laws. 

In this Letter we report a three-dimensional (3D) imag- 
ing study of the microscopic relaxation in a colloidal glass 
under steady shear. The relaxation is nearly isotropic but 
different from that of unsheared colloidal fluids. The in- 
verse relaxation time r" 1 and the diffusion constant D 
show marked shear thinning as a function of the local 
shear rate 7: r" 1 oc D oc 7 8 . We find that this local 
behavior contrasts significantly with the global rheology, 
which shows Herschel-Bulkley behavior. 

We used sterically-stabilized polymethylmethacrylate 
(PMMA) particles (radius a — 850nm, measured by light 
scattering, polydispersity < 10% [17} ) fluorescently la- 
belled with nitrobenzoxadiazole and suspended in a mix- 
ture of cycloheptyl bromide and decalin (viscosity 2.6 



mPa-s) for density and refractive index matching. In 
this medium particles acquire a small charge [18| which 
is largely screened by adding 4 mM tetrabutylammonium 
chloride, giving nearly hard-sphere (HS) behavior, with 
a glass transition at volume fraction 4> g ~ 0.58 (deter- 
mined from mean-squared displacements) (l9j j: we work 
at <f) ~ 0.62 (measured by imaging). The reduced shear 
rate, or Peclet number, is Pe = 4:a 2 j/D = 247Tb, with 
.Do the bare diffusion coefficient and tb — a 2 /6Do= 
1.24 s the Brownian time in our system. Our experi- 
ments cover the range 0.005 < Pe < 1. 

We employ a linear parallel-plate shear cell with a plate 
separation ~ 400 — 800 fim, parallel to ±5^m over a 
~ 200 mm 2 drop of colloid confined between the plates 
by surface tension. We define x, y and z as the velocity, 
vorticity (or neutral) and gradient directions respectively. 
The top plate is driven at 0.05 — 10 /im/s by a mechan- 
ical actuator with magnetic encoder, and steady shear is 
applied up to a total accumulated strain of A7 ~ 1000%. 
Wall slip and wall-induced ordering were prevented by 
coating the slides with 1-3 disordered layers of particles. 
A solvent bath minimized evaporation. 

A 30 x 30 x 15 /an 3 volume in the drop (containing 
N ~ 3000 particles) was imaged from below as a stack of 
75 slices using a fast confocal scanner (VT-Eye, Visitech 
International) and a Nikkon TE Eclipse 300 inverted mi- 
croscope. The scanning of each 3D stack took 1.7 s. Par- 
ticles were located with resolution 5x, 5y ~ 30 nm and 
Sz ~ 90 nm [2p| . Tracking from frame to frame was 
achieved by first subtracting from the raw coordinates 
a time (t) dependent x-displacement profile Ax(z,t), 
measured via correlation analysis of raw images, and 
adding this back after particle tracking. The resulting x- 
displacements over a given time interval dt, {Axi(zi, dt)} 
(i = 1 to N), always have an average linear dependence 
on z. From this we checked that the sample in our im- 
aged volume was indeed subjected to uniform shear, and 
measured the actual (local) shear rate 7, which may differ 
from the applied (global) rate 7 a due to shear localiza- 
tion and the presence of jammed regions. We will return 
to this point; for now we focus on steady states with a 
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FIG. 1: Colloid trajectories for 7 = 9.3 x 10~ 4 s -1 . (a) 
1.5 /im thick slice in the x, z plane for 160 s; The start of each 
trajectory is shown by o, the end by □. The big arrow marks 
the shear direction, (b) As in (a) but in the de-sheared, x,z, 
reference frame, with Xi = Xi — 7 J* Zi(t')dt' . (c) y,z plane 
over 160 s. (d) x,y plane over 160 s. (e) Single trajectory 
in the x, y plane over 800 s. Dotted circles mark rattling 
in several cages (not the particle size), grey dots show the 
locations at t = 0, 200, 400, 600, 800 s. 



linear velocity profile in a region from 15 — 30 /xm above 
the cover slide. When present, strong decay in the shear 
rate occurs at least Az ~ 20a away from imaged regions. 
We also checked, via bond-order analysis [2l|, that shear- 
induced crystallization was absent for our range of 7 [22| • 

Figure [TJa) shows the trajectories in an x, z slice at 
7 = 0.93 x 10~ 3 s _1 . The displacement gradient due 
to shear is evident. To highlight the shear- induced dy- 
namics, we show in Fig. []Jb) the non-affine component 
of the motion obtained by subtracting the uniform shear 
via it — Xi — 7 J* Zi(t')dt' . Considerable shear-induced 
non-affine displacements are seen in this plane as well 
as in the other planes, Figs. [TJc,d). On the time scale 
considered here, these rearrangements are heterogeneous, 
somewhat similar to observations in quiescent concen- 
trated colloidal fluids for cp < <f) g (23[. Zooming in on a 
single particle, Fig. QTe), we observe that its dynamics 
under shear consists of intervals of cage 'rattling', inter- 
rupted by shear-induced plastic cage-breaking events. 

Next, we study the relaxation via the incoherent scat- 
tering function, F s (Q,t) = (cos(Q[yi(t +t)-yi(to)}))i,t , 
at a scattering vector Q = Q m ~ 3.8a -1 where the data's 
structure factor S(Q) shows a peak. In Fig. [2] we show 
selected results for Q \\ y, but the results (not shown) for 
Q || z and x, using the non-affine displacements Xi for 
the latter, are similar. F s for the quiescent glass (7 = 0) 
hardly decays over our observation window, reflecting the 
caging of particles by their neighbors; at longer times we 
observed aging, as in other studies [H HI- The short 
time decay due to initial cage exploration (t < tb [HI, 



FIG. 2: Incoherent scattering functions F s (Q m ,t), with 7 in- 
creasing from right to left. Lines for 7 = 0.93 x 10~ 3 s -1 show 
two curves used in the average with start times io spaced by 
180s. The dashed line schematizes initial relaxation. Inset: 
data collapse using f a (Qm,t/T a ). Line: f s oc exp(— t/r a ). 



dashed line in Fig. [J) is inaccessible to us. At small 7, 
F s at short times still exhibits a plateau, in agreement 
with the caging in Fig.QJc). As 7 increases, this plateau 
shrinks and for the highest 7 it vanishes and likely merges 
with the short time decay. At longer times, F s decays 
strongly for all 7 7^ 0, marking shear-induced structural 
relaxation and cage rearrangements. The structural (a- 
)relaxation time T a , defined by F s (Q m ,t = T a ) = e , 
decreases on increasing 7. Importantly, F s is independent 
of the starting time to (see data for 7 = 0.93 x 10 -3 s _1 ), 
i.e., a stationary state is achieved. 

Our data confirms the theoretically-predicted 'time- 
shear superposition principle' 0, when time is 
scaled by r a , the a-relaxation follows a master curve 
f s (Q,t/T a ), Fig. O inset. As in LJ simulations @, our 
f s is a pure exponential. This differentiates a shear- 
melted glass from a dense HS colloidal fluid at <fr < <f) g 
and 7 = 0, where we find stretched exponential behavior 
for F S (Q > Q m /2), as can also be deduced from [25| . 

Figure [2£a) shows the dependence of T a on the shear 
rate. It exhibits a power law r Q oc 7 _l/ with v = 
0.80 ±0.01 [1^|, independent of the criterion or Q used to 
determine r a . This behavior means that the accumulated 
strain at r a is not constant but varies as jT a oc 7 2 . The 
data are consistent with a schematic model [H for driven 
glasses and also match the 'creep' behavior of a driven 
particle in a correlated random potential 1271. We note 
that an 'cntropic barrier hopping' model [28j, without 
any 'ideal' glass divergences, shows a very similar depen- 
dence of the 'hopping' time on 7. Below we discuss the 
rhcological implications of this behavior. 

Turning to the mean squared displacement (MSD) 
(dy 2 (t)), Fig.[3ja) inset, we see that it exhibits a crossover 
from caged to diffusive motion for -J (dy 2 )/a ~ 0.15 
((dy 2 ) ~ 0.017 £tm 2 ), in reasonable agreement with the 
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FIG. 3: (a) Structural relaxation time (•) and the characteris- 
tic time T2 for the crossover from caged to diffusive behavior 
(■) vs. 7; dashed line: r a oc 7 -0 ' 8 ; full line: T2 oc -y" - 65 . 
Inset: mean square displacement in the vorticity direction 
for shear rates as in Fig. [5] Line: {dy 2 (t)) = 2D y t for 
7 = 1.48 x 10~ 3 s -1 . (b) (•) Diffusion constant D y vs. 
7. Line: D y cx 7 ' 8 . (c) The scaled diffusion constant 
vs. 7 for j = 2;, y, z. Line: the value 
(dy 2 (r a )) — 2/Qm expected for gaussian behavior. 



'Lindemann parameter' measuring the cage rattling at 
the quiescent glass transition [T^]. The long time diffu- 
sion constant D y , Fig. \5j[b), follows the relaxation rate 
D y oc r" 1 oc 7 8 , and woi the shear rate 7. To show 
this more clearly and also address the anisotropy in the 
dynamics, we plot in Fig. [UJc) the product 2DjT a for 
the three directions (j = x, y, z) along with the value 
(dy 2 (T a )) = 2/Q 2 n expected from a gaussian approxi- 
mation F s (Q m ,t) ~ e-Om^WJ/a The value for 
2D y T a agrees well with 2Q^ and this gaussian long 
time behavior also occurs in the other directions [3fj| . 
We again stress the difference with quiescent fluids at 
(f> < 4>g, which always show Dr a < Q~ 2 for Q > Q m . 
Figure [3{c) also shows that the diffusion constants ex- 
hibit only a mild anisotropy: while D z > D x ^ y , the dif- 
ference is < 20%. Similar or even smaller anisotropy has 
been observed in simulations of sheared, glassy systems 
d, Oil and colloidal fluids [3l[ . Isotropic shear- induced 
diffusion is also seen in dilute suspensions [32I ]. However, 
sheared non-Brownian suspensions (Pe — ► 00) show a 
marked anisotropy (D^/D^ z ~ 8) with D°° oc 7. 

As a last characterization of the microscopic dynam- 
ics we study the probability distribution of the displace- 
ments P(dy(t)) and the non-gaussian parameter 0:2,3/ = 
(dy A (t)) /3(dy 2 (t)) 2 — 1. The latter characterizes broad, 
non-gaussian, tails to P(dy(t)), reflecting cage rearrang- 




FIG. 4: (a) Nongaussian parameter 02 (t) of the probability 
distribution P[dy(t)], for several 7. (b) P[dy(t — T2( A /))] for 
the corresponding 7, showing a near collapse of the data (each 
involving > 10 5 displacements). Line: best gaussian fit. 

ing motions such as in Fig. 1(e). Figure[4ja) shows a2, y (i) 
for various 7. It exhibits a peak for t = Ti corresponding 
to the crossover from caged to diffusive behavior in the 
MSD (inset, Fig. EJa)), and vanishes for t > r a . A non- 
zero cxi also suggests cooperative motion, consistent with 
the heterogeneous trajectories for t < r Q in Fig. Q] (b)- 
(d). The peak time follows T2 oc 7 65 , Fig. El^a), some- 
what different from the r Q scaling. More interestingly, 
the distributions P[dy(t = T2(7))] show a near collapse 
for different rates, (Fig. H^b)), despite a slight decrease 
of 0:2(72) with 7. In quiescent systems at (f> < <p g , such 
(near) collapse of P[dy(T2(cf>))] at different <f> is not ex- 
pected since there ui (T2) grows strongly with <f> while 
the MSD at t-i decreases rapidly [23| . 

We now return to the 7 dependence of r a . There 
is currently no firm theoretical basis for relating t q to 
flow properties. Nevertheless, r a is often taken (with 
some simulational evidence [23) as proportional to vis- 
cosity [||, giving an effective stress a = Got q 7 with Go 
a modulus. The resulting 'microscopic' flow curve shows 
a oc 7 2 , Fig. [5] Recent theories 0, Q have argued 
for the existence of a dynamic yield stress at 7 — > + in 
uniform shear. However, our results show no sign of a 
plateau in a for reduced rates down to Pe ~ 0.005. 

Figure [5] shows the experimental global flow curve mea- 
sured with a stress controlled rheometer (AR2000, TA 
Instruments) in cone-plate geometry (diameter 40 mm, 
angle 1°, both surfaces coated with particles). The stress 
a is related to the average shear rate j a by er(7 Q ) = 

+ A'jJ 1 with a dynamic yield stress = 1.36 Pa 
and n = 0.56, similar to previous HS measurements [34[. 

To compare with the microscopic behavior a = GoT a 7, 
in Fig[S] we have chosen Go to optimize agreement be- 
tween a and c(7 a ) at high 7. Clearly, the microscopic 
and macroscopic data disagree. Some discrepancy may 
be due to the fact that, for Pe > 1, hydrodynamic effects 
render the relation a oc Taj less valid. More importantly, 
discrepancy could arise from the presence of shear local- 
ization, e.g. due to the existence of a static yield stress 
0, [l6[ ■ We have already mentioned that in our parallel- 
plate shear cell, a global shear rate j a typically corre- 
sponds to a jammed region (7 ~ 0) coexisting with a flow- 
ing region with 7 > j a . Preliminary flow imaging inside 
our rheometer shows that shear localization also occurs 
in the cone-plate geometry, and sets in for j a < 10 -2 s _1 
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FIG. 5: (□) Local 'flow curve' a = Gojr a , with Go = 8.5 Pa, 
vs. 7 (dashed line: a oc 7°' 2 ) compared with the macroscopic 
flow curve cr(7a) measured in cone-plate geometry (•). Full 
line: fit to the Herschel-Bulkley model a = 1.36 Pa + A^°- 56 . 

[35j ] . In LJ simulations differences between 7 and 
7o could indeed explain small deviations between local 
and global rheology. But the rather larger deviations in 
Fig. O cannot be explained by this argument. Possibly, 
the relation a = Got q 7 is an oversimplification [3(| and 
instead we may need to invoke analogies with 'force chain' 



dominated systems to make progress. Indeed, our global 
shear profile 7(2), which exhibits a smooth rather than 
a step-like decay of 7(2) to zero (data not shown), has 
similarities with velocity profiles in granular matter [37| . 

Concluding, we have studied the 3D particle dynamics 
in a HS colloidal glass under steady shear by fast confo- 
cal imaging. Shear occurs in 'fluidized' bands where col- 
loids show nearly isotropic 'cage breaking' and exponen- 
tial relaxation, in contrast to the stretched-exponential 
dynamics in dense colloidal fluids. The relaxation rate 
scales as a power of the local shear rate: r" 1 oc 7 8 . 
The 'naive' microscopic flow curve deduced from this re- 
sult differs from the global, Herschel-Bulkley, rheology. 
These and other recent results HH show the potential 
of fast 3D imaging to address fundamental questions in 
non-equilibrium physics. 
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